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HP-41CX  PROGRAMS 
FOR 

HgCdTe  DETECTORS  AND  IR  SYSTEMS 


I.  INTRODUCTION 

Programs  have  been  written  for  the  HP-41  CX  which  aid  in  the  analysis  of  HgCdTe  detectors,  focal 
planes,  and  infrared  systems.  They  have  been  installed  as  a  “bundle”  in  an  HP-41CX  and  a  “WALL” 

(write-all)  operation  has  been  performed,  saving  the  group  of  programs  and  machine  status  to  magnetic 
cards.  They  have  been  written  to  run  in  a  basic  HP-41  CV  or  HP-41  CX  with  no  card  reader  or  additional 
ROMs  required.  An  HP-41C  with  less  than  maximum  RAM  will  be  unable  to  read  this  WALL 
configuration  and  will  be  unable  to  hold  all  the  programs  in  memory  at  the  same  time,  but  a  subset  can 
be  loaded  by  hand  from  the  enclosed  listings. 

The  purpose  of  this  report  is  to  describe  the  status  of  the  HP-41CX  after  the  WALL  cards  have  been 
read,  and  to  provide  the  user  with  information  needed  to  use  the  programs. 

It  is  assumed  that  the  user  has  some  familiarity  with  the  operation  and  programming  of  the  HP-41 
family  of  pocket  calculators.  However,  the  description  will  begin  with  a  brief  review  of  the  more  pertinent 
characteristics  and  procedures. 

The  HP-41CX  can  have  many  different  unique  programs  in  memory  at  the  same  time.  Jumps  to 
perform  a  subroutine  or  to  alter  program  flow  are  referred  to  labels  (LBLs)  which  identify  the  target 
location.  There  are  a  limited  number  of  short  (quick)  labels  available,  and  when  writing  a  program,  it’s 
important  to  keep  track  of  label  usage.  The  HP-41CX  provides  the  capability,  when  desired,  of  separating 
programs  with  an  END  statement.  The  short  labels  cannot  be  reached  by  GTO  (go  to,  jump)  or  XEQ 
(execute,  gosub)  commands  beyond  an  END  statement.  This  means  that  such  labels  are  local  labels  and 
don’t  conflict  with  identical  labels  in  other  programs  separated  by  an  END  statement.  There  are, 
however,  global  labels  available  which  permit  the  user  to  go  directly  to  a  program  from  anywhere  in 
RAM.  These  also  permit  one  program  to  GTO  or  XEQ  another  program,  even  though  they  are  separated 
by  an  END  statement.  Performing  a  CATALOG  1  causes  the  HP-41  CX  to  list  the  global  labels.  The 
listing  may  be  paused  with  the  R/S  key,  and  you  may  step  forward  or  backward  through  the  listing  with 
the  SST  and  BST  keys.  You  may  put  yourself  in  a  desired  program  by  listing  CATALOG  1,  pausing  with 
the  desired  program  in  the  display,  and  pressing  the  ♦-  key. 

The  global  labels  of  the  programs  included  in  this  WALL  configuration  are  listed  on  the  next  page. 

Global  labels  need  not  necessarily  be  separated  by  END  statements,  and  some  of  the  labels  in  this  group 
of  programs,  in  fact,  are  not.  Indented  programs  are  included  in  the  same  local  area  with  program  above. 

LABEL  DESCRIPTION 

FLAG  Puts  calculator  flags  back  to  a  default  status. 

NEdt  Noise  equivalent  delta  temperature  of  an  IR  system,  and  photons  collected  in  an  integration  time. 

DET  RNG  Estimates  detection  range  of  a  target. 

D*  D-Star  (PV,  Johnson,  and  BLIP),  RoA  BLIP,  D*lambda  to  D*BB  conversion  factor. 

QBB  Blackbody  photon  flux. 

BB  General  purpose  blackbody  calculations. 

SV  General  purpose  root  finder  used  by  other  programs. 

GR-DIFF  RoA,  temperature,  cutoff  interpolator /extrapolator. 

HCT  X  of  Hg(x)Cd(l-x)Te  from  cutoff  and  temperature. 

EG  Bandgap  from  x  and  temperature. 

NI  Intrinsic  carrier  concentration  from  material  parameters. 

R0A  RoA  product  from  material  parameters. 

CONSTNT  Useful  physical  constants. 

INJEFF  Direct  Injection  injection  efficiency. 


The  following  USER  key  definitions  are  also  included.  (In  USER  mode,  the  following  keys  perform  the 
indicated  operation.) 


PRESSING 

WILL  RESULT  IN 

TAN 

PACK 

SHIFT  TAN 

CLRG  (clear  all  data  registers) 

SHIFT  R/S 

Run  the  FLAG  program 

Twenty-six  registers  have  been  reserved  for  data  storage  (registers  00  through  25.)  Thirty- five  registers 
of  the  319  available  in  an  empty  machine  are  unused  and  available  for  either  additional  programming  or 
additional  data  storage.  Each  register  has  room  for  approximately  five  program  steps.  Should  more  room 
be  needed  for  a  different,  lengthy  program,  and  if  you  are  using  an  HP-41CX  or  an  HP-4lC(V)  with  an 
extended  memory  module,  you  may  SAVEP  any  program  in  RAM  to  EXTENDED  MEMORY  (which  is 
similar  to  a  RAM-DISK)  and  then  clear  the  program  from  RAM.  This  may  be  repeated  as  necessary  with 
other  programs  until  EXTENDED  MEMORY  is  full.  Programs  in  EXTENDED  MEMORY  cannot  be 
run,  but  you  may  GETP  them  back  into  RAM,  if  there  is  enough  RAM  available. 

Of  the  26  data  registers  available,  only  the  first  15  STO  (store)  and  RCL  (recall)  statements  are  quick 
access  (one-byte  statements  instead  of  two-byte  statements.)  You  can  ignore  this  fact,  but  it  is  the  reason 
that  the  first  15  registers  are  used  most  often.  Since  there  are  more  than  15  different  physical  parameters 
involved  in  the  above  programs  and  some  data  registers  are  needed  for  bookkeeping  or  temporary  scratch 
storage,  there  is  not  a  firm  relationship  from  global  program  to  global  program  between  any  physical 
parameter  and  any  data  register.  You  can  rest  assured  that  the  right  data  are  in  the  right  place  at  the  right 
time  when  running  a  given  program.  BUT  don’t  assume  that  data  entered  while  running  one  program  will 
be  correctly  placed,  should  you  go  to  another  global  program  which  requires  the  same  data.  It  is 
advisable  to  perform  a  CLRG  (SHIFT-TAN  in  USER  mode)  and  reenter  all  necessary  data  when  moving 
to  a  new  program. 


There  is  a  similar  caution  concerning  the  HP-41CX  flags.  There  are  56  flags  which  are  used  by  the 
HP-41CX,  30  of  which  are  user  alterable.  The  status  of  the  first  five  are  indicated  on  the  LCD  display. 
Since  the  calculator  is  slow  (as  computers  go),  some  calculations  are  done  once,  the  data  are  stored,  and  a 
flag  is  set  to  let  the  program  know  that  the  calculation  need  not  be  done  again,  unless  certain  parameters 
are  altered.  If  these  parameters  are  altered,  the  flag  is  automatically  cleared  so  that  a  new  calculation  will 
be  performed  when  required.  Flags  are  used  to  control  program  flow  for  other  reasons  as  well.  To  keep 
the  user  informed,  only  the  first  five  flags  (which  are  visible  in  the  display)  are  used  in  this  fashion. 
Different  global  programs  may  use  the  same  flag  to  mean  different  things,  so  it  is  advisable  when  moving 
to  a  new  global  program  (or  just  starting  a  program)  to  check  to  see  that  the  five  displayed  flags  (0 
through  4)  are  clear.  If  you  can’t  see  them,  they’re  clear.  Running  the  FLAG  program  (SHIFT-R/S  in 
USER  mode)  will  insure  that  all  flags  are  clear.  This  will  leave  you  in  the  FLAG  program,  and  then  you 
will  have  to  go  to  the  program  you  wish  to  run. 

Running  a  program  generally  proceeds  as  follows.  With  the  exception  of  FLAG  which  may  be  run  at 
any  time  by  pressing  SHIFT  R/S,  the  user  must,  first,  position  himself  within  the  program  of  interest. 
This  may  be  done,  as  described  above,  by  using  CATALOG  1  and  pressing  R/S  to  stop  the  catalog  scan 
with  the  desired  program  name  in  the  display,  or  by  going  directly  to  the  program  of  interest  with  the 
GTO  command.  For  example,  to  run  the  blackbody  program  (BB),  press  keys  as  follows:  SHIFT  GTO 
ALPHA  BB  ALPHA. 


Make  sure  that  flags  0  through  4  are  clear.  It  may  be  wise  to  do  a  CLRG  to  increase  the  likelihood 
that,  should  you  fail  to  provide  all  the  necessary  data  to  the  program,  an  error  will  occur  or  invalid 
results  will  be  clearly  invalid. 

The  programs  generally  require  data  to  be  input  before  any  results  may  be  computed.  The  top  row  of 
keys  (A,  B,  C,  D,  E)  and  their  SHIFTed  counterparts  (a,  b,  c,  d,  e)  are  used  to  communicate  with  a 
program.  This  is  because,  WHEN  IN  USER  MODE,  pressing  one  of  these  keys  will  begin  program 
execution  at  the  label  (LBL)  which  matches  the  key.  These  labels  are  local.  If  no  such  label  exists  in  the 
current  program  area,  then  the  function  identified  on  the  keyboard  is  executed  instead.  For  example, 
suppose  the  current  program  area  contains  the  instruction  sequence: 


LBL  d 
5 


With  USER  mode  active,  pressing  SHIFT  D  will  cause  the  program  to  being  running  at  LBL  d,  where  it 
will  multiply  the  contents  of  the  X  register  (the  display)  by  S,  then  add  1  and,  finally,  stop  showing  the 
answer  again  in  the  X  register. 

The  programs  are  written  so  that  pressing  certain  keys  will  cause  data  which  have  been  entered  into  the 
X  (or  X  and  Y)  (or  X,  Y,  and  Z)  register  to  be  stored  away  for  later  use.  Pressing  certain  keys  (perhaps 
after  some  additional  data  entry)  will  cause  desired  calculations  to  be  performed  and  the  results  displayed. 
Which  keys  do  what  varies  from  program  to  program,  and  the  functions  are  not  indicated  on  the 
keyboard.  You  must,  therefore,  either  remember  the  key  definitions  or  refer  to  a  document  such  as  this 
report.  It  is  common,  when  an  individual,  dedicated  program  is  stored  to  a  magnetic  card,  to  write  the 
definitions  for  each  key  on  that  magnetic  card.  The  card  is  conceptually  divided  into  five  parts 
horizontally  representing  the  five  keys  in  the  top  row,  and  two  parts  vertically  representing  the  unshifted  and 
shifted  labels.  The  upper  half  represents  the  shifted  keys.  Although  this  has  not  been  done  with  the 
WALL  cards,  the  same  convention  will  be  followed  by  including  in  each  program  description  a  drawing 
based  on  this  layout  which  defines  the  operation  of  each  key. 

It  is  somewhat  inconvenient  to  use  a  single  key  to  enter  multiple  data,  but  it  is  necessary,  since  there  is 
often  more  data  required  for  a  calculation  than  there  are  keys.  Referring  to  one  of  these  drawings, 

ENTER  means  press  the  enter  key  on  the  keyboard  after  keying  in  the  requested  data.  This  pushes  the 
data  just  entered  up  into  the  Y  register  and  permits  the  user  to  enter  additional  data  in  the  X  register. 

Another  ENTER  pushes  the  contents  of  Y  into  Z,  X  into  Y,  so  additional  data  may  be  placed  in  X.  After 
keying  in  all  the  requested  data  for  a  key,  press  that  key  so  that  the  program  will  store  it.  Continue  from 
key  to  key,  entering  in  all  data  required  to  compute  the  parameter  desired.  Note  that  some  programs 
calculate  more  than  one  parameter  and  that  some  parameters  may  not  require  that  all  the  data  indicated 
on  the  drawing  be  input.  Which  parameters  require  which  data  comes  from  experience  or  referral  to  this 
report. 

Once  all  the  necessary  data  have  been  keyed  and  stored,  pressing  the  FIND  key  for  the  desired 
parameter  will  initiate  the  calculation  to  produce  the  answer  for  that  parameter.  Obviously,  pressing 
FIND  before  all  the  necessary  support  data  have  been  entered  will,  very  likely,  produce  invalid  results.  It 
may  even  produce  an  operational  error;  e.g.,  divide  by  zero,  causing  the  calculation  to  terminate.  Under 


most  circumstances  no  harm  is  done  when  this  occurs.  However,  NOTE  that  some  programs  call  other 
programs.  If  the  error  occurs  in  the  called  program,  you  will  be  left  in  that  program,  not  in  the  program 
from  which  you  started.  It  is  advisable,  should  such  an  improper  termination  occur,  to  re-position 
yourself  via  CATALOG  1  or  GTO  to  the  program  of  interest.  It  is  probably,  also,  wise  to  clear  any  set 
flags. 

If  incorrect  data  have  been  keyed  and  stored,  it  may  be  corrected  by  performing  the  correct  data  entry 
sequence  for  the  key  dedicated  to  that  data.  This  means  that  you  will  have  to  reenter  all  the  data 
associated  with  that  key.  This  is  also  true  if  you  wish  to  change  any  data,  either  before  or  after  FINDing 
a  result.  To  recalculate  a  result  with  different  initial  conditions,  you  need  only  reenter  the  data  associated 
with  the  key  whose  data  you  wish  to  alter.  Then,  press  the  appropriate  FIND  key  for  the  result  desired. 

What  follows  is  user  documentation  and  description  for  each  of  the  programs. 


II.  PROGRAM  DOCUMENTATION 


1.  FLAG.  This  is  a  very  short  program  designed  specifically  to  restore  all  55  flags  of  the  HP-41CX  to 
default  status.  It  dears  Hags  0  through  25,  sets  the  display  to  FIX  2  format,  sets  DEG  mode,  calls  for  a 
period  as  the  decimal  point,  and  calls  for  commas  between  digit  groupings  of  three.  It  alters  no  data 
registers.  Ti  is  program  may  be  quickly  run  by  pressing  SHIFT  R/S  in  the  USER  mode  if  FLAG  has  been 
assigned  to  this  key. 

The  FLAG  program  on  the  WALL  cards  contains  undocumented  program  steps  and  should  not  be 
altered  or  single-stepped  in  RUN  mode. 

Also,  included  is  a  listing  for  FLAG  which  simulates  the  short  version  and  which  contains  no 
undocumented  operations.  This  version  is  significantly  slower  than  the  short  version,  but  it  may  be  keyed 


by  hand  and  needs  no  precautions. 


PROGRAM  LISTINGS  FOR  FLAG 


Two  listings  are  given.  The  first  is  very  fast  and  short  but  contains  undocumented  code  which  can  only 
be  created  by  procedures  unavailable  to  most  users.  It  is  included  on  the  WALL  cards.  The  second  is 
longer  and  significantly  slower  but  may  be  easily  keyed  by  hand. 


81/10/1986 


elm  -aflc- 

02  *<x»* 

83  fiSTO  d 

84  CF  83 

85  . END. 


01*LBL  -FIRG- 
82  .025 
03»LBL  90 
04  CF  IND  X 
05  ISG  X 
06  GT0  06 
07  SF  26 
08  SF  27 
09  SF  28 

10  SF  29 

11  FIX  2 

12  DEG 

13  RTH 

14  END 


2.  NEdT,  DET  RNG. 


Lambda  Initial 
ENTER 
Lambda  Final 
(mM) 

SHIFT  R 

.  /#  cold  shld 
ENTER 

F/#  signal 

SHIFT  B 

Bkgnd  temp 
ENTER 
Detector  temp 
(Kelvin) 

SHIFT  C 

lnteg.  time 
(sec) 
ENTER 
Quant,  eff. 
(percent) 

SHIFT  D 

Detector  area 
(sq.  cm.) 

SHIFT  E 

FIND 

FIND 

Target  dT  (K) 

S/N  criteria 

RoA 

Photons  per 

ENTER 

ENTER 

(ohm-cm  2) 

NEdT 

integration 

Target  Area 

Atm  transmissn 

(K) 

time. 

(sq.  m) 

(percent/km) 

ENTER 

FIND 

Focal  length 

Detection  Rng 

(mm) 

(km) 

u  n  _ 

\ 

-1 

Below  is  a  table  showing  the  data  required,  registers  used,  sample  test  data,  and  results. 


REQ’D  FOR  REQ’D  FOR  REQ’D  FOR  REGISTER  TEST 
NEdT  PHOTONS  DET  RNG  USED  DATA 


LAMBDA  INITIAL 

X 

X 

X 

1 

3 

LAMBDA  FINAL 

X 

X 

X 

5 

F/#  COLD  SHIELD 

X 

X 

17* 

1 

F/#  SIGNAL 

X 

X 

X 

4 

1.4 

BACKGROUND  TEMPERATURE 

X 

X 

X 

3 

DETECTOR  TEMPERATURE 

X 

X 

14 

INTEGRATION  TIME 

X 

X 

X 

5 

4E-5 

QUANTUM  EFFICIENCY 

X 

X 

X 

18 

DETECTOR  AREA 

X 

X 

X 

2 

IE-5 

RoA 

X 

X 

15 

TARGET  DELTA  TEMP 

X 

16 

3 

TARGET  AREA 

X 

11 

3 

SYSTEM  FOCAL  LENGTH 

X 

12 

S/N  CRITERIA 

X 

13** 

ATM.  TRANSMISSION  (PERCENT/KM) 

X 

20*** 

70.71 

RESULTING  NEdT  =0.1 551  IK  Other  registeres  used: 

RESULTING  SIGNAL  PHOTONS  =  4.02  E5  0,6, 7, 8,9 

RESULTING  DETECTION  RANGE  =  3.885  km 

*  This  register  actually  contains  the  ratio  of  the  f/numbers  squared. 

••  This  register  used  for  scratch  as  well,  so  S/N  not  saved. 

*•*  SQR  (transmission  percent/ 100)  actually  stored. 


This  program  finds  (1)  Noise  equivalen,  delta  temperature  (NEdT),  (2)  number  of  photons  collected  by 
a  detector  (or  an  equivalent  detector,  for  TDI)  through  the  system  in  an  integration  time,  and  (3) 
detection  range  for  a  target. 

The  NEdT  for  an  IR  system  is  calculated  from  certain  system  and  detector  parameters.  It  is  a  relatively 
abbreviated  calculation,  so  the  results  are  only  approximate.  Assumptions  made  include  that  the  noise 
comes  only  from  the  detector  or  the  background,  that  the  noise  bandwidth  is  the  reciprocal  of  the 
integration  time,  that  the  quantum  efficiency  is  constant  with  respect  to  wavelength,  and  that  the  detector 
noise  is  due  to  shot  noise  at  zero  bias. 

The  calculation  performed  is: 

NEdT  =  SQR(2*PC  +  (4kTt/Ro/eA2)), 

SC 

where  Q.E.  is  the  detector  quantum  efficiency; 
t  =  detector  integration  time; 

T  =  detector  temperature; 

Ro  =  detector  zero  bias  dynamic  resistance; 

A  =  detector  active  area;  and 

PC  is  the  charge  absorbed  by  the  detector/t  due  to  photons. 

or,  PC  =  Q  A  (Q.E.)  t 

(2  *  (cold  shield  F/0))A2 

where  Q  =  photon  flux  rate  from  cuton  to  cutoff/cmVs,  and 

SC  is  the  signal  charge  absorbed  by  the  detector/t  per  degree  change  in  scene  temperature. 

or,  SC  =  dQ  A  (Q.E.)  t  , 

dT  (2  *  (signal  F/#))A2 

where  dQ  =  change  in  signal  flux  per  change  in  scene  temp. 
dT 

The  photons  per  integration  time,  found  by  pressing  key  B,  is  the  number  of  SIGNAL  photons 
collected  through  the  optics  of  the  indicated  system  (using  the  signal  T/N)  which  produces  signal  carriers. 
For  TDI  systems  the  integration  time  entered  should  be  the  integration  time  per  detector  times  the  number 
of  detectors  in  TDI,  so  that  the  number  of  carriers  accumulated  is  the  number  for  the  TDI  column.  To 
find  the  number  of  background  carriers  collected,  set  the  signal  F/0  equal  to  the  background  F/0. 

The  detection  range  is  calculated  by  finding  the  range  at  which  the  image  of  the  target  is  sufficiently 
smaller  than  the  detector,  so  that  the  S/N  is  as  desired.  The  signal  is  attenuated,  based  upon  the  entered 
atmospheric  transmission  per  kilometer.  The  SV  (solve)  program  is  called  to  perform  the  iteration  to  a 
solution. 

S/N  =  (transmissionf(range)  x  (target  delta T)  x  (target  area)  x  fA2 

NEdT  A  Ranged 


PROGRAM  LISTING  LOR  NFdT,  DET  RNG 


01/18/1986 

51  RCL  13 

52  * 

53  RCL  97 

01«LBL  ’BET  RNG* 

54  - 

02*L8L  C 

55  RTN 

03  STO  12 

56*LBL  *NEdT 

04  RON 

57*LBL  A 

05  STO  11 

58  .5 

06  RON 

59  ST*  93 

07  STO  16 

60  SF  00 

08  RTN 

61  XEQ  38 

09*L8L  D 

62  STO  96 

10  100 

63  1 

11  / 

64  ST-  93 

12  SORT 

65  XEQ  08 

13  STO  20 

66  CF  30 

14  XOY 

67  .5 

15  STO  13 

63  ST+  03 

16  50  *NEdT* 

69  RDN 

17  RCL  13 

70  RCL  06 

18  * 

71  XOY 

19  RCL  16 

72  ST-  06 

20  XOY 

73  * 

21  *N0  DETECT* 

74  RCL  17 

22  X>Y? 

75  / 

23  PROMPT 

76  215  E13 

24  / 

77  RCL  14 

25  Ra  12 

78  * 

26  Xt2 

79  RCL  15 

27  RCL  11 

38  / 

28  * 

81  RCL  02 

29  RCL  02 

82  * 

38  / 

83  RCL  05 

31  * 

84  * 

32  SORT 

85  + 

33  1  E4 

86  SORT 

34  / 

87  RCL  06 

35  STO  13 

83  / 

36  SF  10 

89  FIX  5 

37  *R* 

99  *NEdT=* 

38  ASTO  06 

91  ARCL  X 

39  8 

92  flVIEH 

40  XOY 

93  RTN 

41  FIX  3 

94*LBL  a 

42  XEQ  *SY* 

95  STO  19 

43  CLA 

96  XOY 

44  ARCL  X 

97  STO  01 

45  *H(H* 

93  RTN 

46  PROMPT 

99*L6L  b 

47»L8L  *R* 

188  x< r 

48  RCL  20 

101  GTO  99 

49  XOY 

102  STO  84 

50  YtX 

193  / 

104  Xt2 

157  Xt2 

105  STO  17 

158  / 

186  RTN 

159  + 

107*LBL  c 

160  RCL  08 

108  STO  14 

161  9 

109  XOY 

162  * 

110  STO  03 

163  6 

111  RTN 

164  + 

112*LBL  d 

165  RCL  08 

113  STO  18 

166  * 

114  XOY 

167  2 

115  STO  85 

163  + 

116  RTN 

169  27 

117»LBL  e 

170  / 

118  STO  82 

171  RCL  07 

119  RTN 

172  XT2 

129*LBL  E 

173  / 

121  STO  15 

174  RCL  07 

122  RTN 

175  / 

123*LBL  B 

176  + 

124  CF  88 

177  RCL  03 

125»LBL  88 

173  Xt2 

126  SF  91 

179  * 

127  RCL  01 

188  RCL  83 

123  XEQ  81 

131  * 

129  STO  80 

132  1581  E5 

138  RCL  18 

183  * 

131*LBL  01 

184  RCL  18 

132  14388 

185  * 

133  XOY 

186  RCL  04 

134  / 

187  Xt2 

135  RCL  83 

188  / 

136  / 

189  RCL  82 

137  STO  88 

199  * 

138  2 

191  RCL  05 

139  * 

192  * 

140  RCL  88 

193  FS?C  81 

141  * 

194  RTN 

142  2 

195  RCL  09 

143  + 

196  - 

144  RCL  88 

197  FS?  98 

145  EtX 

198  RTN 

146  STO  07 

199  SCI  2 

[47  / 

208  CLA 

1«8  RCL  88 

201  ARCL  X 

149  XT2 

282  *h  PHOT- 

158  RCL  93 

283  AVI  EH 

151  * 

204  RTN 

152  .5 

265  END 

153  * 

154  2 

155  / 

156  RCL  87 
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3.  D*,  QBB. 


aBHi 

p  Lambda  Initial 

ENTER 
Lambda  Final 
(pM) 

SHIFT  R 

Quant.  Eff. 
(percent) 
ENTER 
Field  of  View 
(full  angle) 
(degrees) 

SHIFT  B 

Bkgnd  temp 
ENTER 
Detector  temp 
(Kelvin) 

SHIFT  C 

RoA 

(ohm-cm  2) 

SHIFT  D 

FIND 

RoA  BLIP 

SHIFT  E 

FIND 

FIND 

FIND 

source  temp 
FIND 

FIND 

photon  flux 

D*  (PV) 

D* (BLIP) 

D*  (JOHNSON) 

D*lambda  to 
D*BB 
conversion 
factor 

with  Q.E. 
and  FOV 

_ _ R _ 

.  B  . 

c 

0 

\ 

1 

UJ 

mkhmhhII 

Below  is  a  table  showing  the  data  required,  registers  used,  sample  test  data,  and  results. 


REQ’D  REQ’D 
FOR  FOR 
•BLIP  D*JHNSN 


REQ’D  REQ’D  REQ’D  RGSTR  TEST 
FOR  FOR  FOR  USED  DATA 
C.F.  P.  FLUX  RoABL 


LAMBDA  INITIAL 
LAMBDA  FINAL 
QUANTUM  EFFICIENCY 
FIELD  OF  VIEW 
BACKGROUND  TEMP. 
DETECTOR  TEMPERATURE 
RoA 

SOURCE  TEMPERATURE 


RESULTING  D*PV 

=  1.93  EJO 

OTHER  REGISTERS  USED: 

RESULTING  D*BLIP 

=  4.82  E10 

0,6,7,8,15 

RESULTING  D* JOHNSON 

=  2. II  E10 

RESULTING  CONVERSION  FACTOR 
RESULTING  PHOTON  FLUX 

=  5.60 

=  1.65  E17/CM*2/SEC 

RESULTING  RoA  BLIP 

=  0.522 

•Source  temperature  is  placed  temporarily  in  register  12;  background  temperature  is  restored  to  register 

12,  at  the  end  of  the  calculation. 

MM 


This  program  contains  routines  to  find  six  different  parameters  associated  with  photovoltaic  detector 
D-star. 

(1)  D* JOHNSON  is  the  detectivity  based  on  the  assumption  that  the  detector  shot  noise  at  zero  bias  is  the  only 
noise  source.  It  is  computed  from: 

D*JOHNSON  =  e  (Q.E.)  (cutoff)  SQR(Ro*A/(4kT)). 

he 

(2)  D*BLIP  is  the  detectivity  based  on  the  assumption  that  the  background  noise  is  the  only  noise  source. 

It  is  computed  from: 

D*BLIP  =  (cutoff)  SQR((Q.E.)/(photon  flux)). 

SQR(2)  h  c  SIN(FOV/2) 

(3)  D*PV  is  the  detectivity  which  includes  both  detector  and  background  noise  sources  and  is  computed 
from  the  above  as: 


D*PV  = 


(D*  JOHNSON)  (D*BLIP) 


SQR((D*JOHNSON)A2  +  (D*BLIP)A2) 


(4)  Blackbody  photon  flux  is  the  photon  collection  rate  per  second  per  square  centimeter  of  detector  area 
for  the  given  quantum  efficiency  and  field  of  view.  It  is  found  by  twice  calling  a  routine  which  calculates 
the  flux  from  lambda  =  0  up  to  lambda.  Lambda  is  first  set  to  cutoff,  then  to  cuton,  and  the  results  are 
subtracted.  The  basic  routine  which  is  called  is: 


Q  =  aSUM(XA2  +  2*X/n  +  2/nA2), 
n  EXP  (nX) 


for  n  =  1 ,2,3  and  X  =  lambda,  and 


where  a  =  2  PI  c  (kT/(hc))A  3. 


(5)  The  conversion  factor  is  the  factor  which  converts  a  blackbody  D*  to  a  lambda  peak  D*.  It  is  a 
theoretical  value  which,  here,  assumes  that  the  quantum  efficiency  of  the  detector  is  constant  between 
cuton  and  cutoff.  For  this  calculation,  it  is  assumed  that  the  peak  response  is  at  the  cutoff.  It  is 
calculated  from  (where  sigma  is  the  Stefan-Boltzmann  constant): 


C.F.  =  (photon  flux)  (b/a),  where  a  is  as  above,  and 


b  = 


2  PI  kA3 


c  (sigma)  (source  temperature)  h  A  2  (lambda  cutoff) 


(6)  RoA  BLIP  is  the  value  of  the  RoA  for  which  D*  JOHNSON  and  D*BLIP  equal  one  another.  It  is, 
therefore,  the  RoA  for  which  the  D*PV  is  .7071  of  D*BLIP.  Requesting  this  calculation  presents  the 
answer  in  the  display  and,  also,  stores  this  value  in  register  9.  This  means  that  further  calculations  for  D* 
will  be  based  on  this  RoA. 
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PROGRAM  LISTING  FOR  D*,  QBB 


01/10/1986 


01H.BL  e 
02  1 

03  STO  09 
04  XEQ  B 
85  XEQ  C 
06  / 

07  Xt2 
08HBL  d 
09  STO  09 

10  RTH 

1 1  ♦LBL  a 

12  CF  00 

13  STO  10 

14  XOY 

15  STO  05 

16  RTH 
17*LBL  b 

18  CF  80 

19  STO  13 

20  XOY 

21  STO  11 

22  RTH 
23HJ5L  c 

24  CF  00 

25  STO  14 

26  XOY 

27  STO  12 

28  RTH 
29HJM.  *«** 
30*LBL  fl 

31  XEO  B 

32  XEQ  C 

33  R-P 

34  STO  08 

35  P-R 

36  * 

37  RCl  08 

38  / 

39  RTH 
48»LBl  E 

41  SCI  2 

42  RCl  15 

43  FS?  09 

44  RTH 

45  RCl  05 

46  XEQ  08 

47  STO  06 

48  RCl  18 


49  XEQ  08 

50  RCL  06 

51  - 

52  FS?  62 

53  RTH 

54  RCl  13 

55  2 

56  / 

57  SIH 

58  Xt2 

59  * 

60  RCl  11 

61  X 

62  STO  15 

63  SF  80 

64  RTH 
654.BL  B 

66  XEQ  E 

67  SORT 

68  RCL  11 

69  XOY 

70  / 

71  RCl  10 

72  * 

73  356  E14 

74  * 

75  RTH 
76tlBL  C 

77  SCI  2 

78  RCl  09 

79  RCl  14 
88  / 

81  SORT 

82  RCl  10 

83  * 

84  RCl  11 

85  * 

86  10854  E5 

87  * 

88  RTH 
89*LBL  -QBB’ 
9041BI  08 

91  3 

92  STO  07 

93  CIX 

94  14388 

95  XOY 

96  / 

97  RCL  12 

98  / 

99  STO  88 
180  CLX 


101*LBL  00 
182  2 

103  RCL  07 

104  / 

105  RCL  08 

106  ♦ 

107  RCl  08 

108  * 

109  2 

110  RCl  87 

111  Xt2 

112  / 

113  + 

114  RCl  07 

115  / 

116  RCL  08 

117  RCl  87 

118  » 

119  EtX 

120  / 

121  + 

122  BSE  07 

123  GTO  80 

124  FS?  02 

125  RTH 

126  RCl  12 

127  Xt2 

128  RCl  12 

129  * 

138  6324  E7 

131  * 

132  ♦ 

133  RTH 
134*LBl  B 

135  XO  12 

136  STO  08 

137  CF  08 

138  SF  92 

139  XEQ  E 
148  CF  02 

141  2216 

142  * 

143  RCl  10 

144  / 

145  RCL  12 

146  / 

147  RCL  80 

148  STO  12 

149  RBH 

150  l/X 

151  EHB 
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4.  BB. 


■'.mKTW'. 


INITIALIZE 


SHIFT  R 

PREFIX 


FIND 

or 

ENTER 
Blackbody 
Temp.  (K) 

SHIFT  B 


PREFIX 


FIND 

or 

ENTER 
Cutoff  O^M) 


SHIFT  C 


PREFIX 


COMPARE 
Lambda  Max 
versus 
BB  Temp, 
(maximized 
photon  flux) 

SHIFT  D 


FIND 


COMPARE 
Lambda  Max 
versus 
BB  Temp, 
(maximized 
power) 

SHIFT  E 

FIND 


d  integral  to  integral  to  Photon  Flux  Power  Flux 

dT  infinity  cutoff 

wrt  lambda 


CHOICE 


Below  is  a  table  showing  sample  test  data  and  results. 

LAMBDA  FINAL  =  5  juM  (not  required  for  integral  to  infinity) 
BLACKBODY  TEMP.  =  300K  REGISTERS  USED:  0  THRU  15 


RESULTS 


PREFIXES  PRESSED 

Q  (photons) 

W  (watts) 

NONE 

2.058  E16 

8.175  E-4 

/cmVs/jiM 

/cm'/jiM 

INTEGRAL  TO  INFINITY  (ITI) 

4.105  E18 

4.592  E-2 

/cm  Vs 

/cm! 

INTEGRAL  TO  CUTOFF  (ITC) 

1.320  E16 

5.901  E-4 

/cmVs 

/cmJ 

DERIVATIVE  WRT  TEMP  (DWT) 

6.580  E14 

2.614  E-5 

/cmVs/uM/K 

/cmVfiM/K 

DWT,  ITI 

4.105  E16 

6.123  E-4 

/cmVs/K 

/cmVK 

DWT,  ITC 

4.749  E14 

2.149  E-5 

/cmVs/K 

/cmVK 

The  COMPARE  functions  permit  you  to  key  in  either  wavelength  or  temperature  and  calculate  the 
other  (temperature  or  wavelength)  for  maximum  flux.  For  example,  keying  in  300  and  pressing  SHIFT  D 
will  produce  the  value  12.23.  This  means  that  for  a  300  K  blackbody,  the  wavelength  of  peak  photon  flux 
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density  is  12.23  *iM.  Or  conversely,  a  blackbody  whose  peak  photon  flux  density  is  at  300  mM  wavelength 
has  a  temperature  of  12.23  K.  The  operation  for  SHIFT  E  is  the  same,  except  the  criterion  is  the 
wavelength  for  peak  power  density.  Keying  300  and  SHIFT  E  produces  the  value  9.659.  NOTE:  The  data 
entered  prior  to  a  COMPARE  operation  does  not  alter  any  values  that  have  been  stored  for  wavelength  or 
temperature  via  SHIFT  B  or  SHIFT  C. 

This  program  uses  the  keys  a  little  differently.  First,  before  doing  anything  else,  you  must  initialize  the 
program  by  pressing  SHIFT  A.  This  configures  the  HP41  data  registers  and  flags  for  the  rest  of  the 
program.  It  also  is  used  to  clear  partial  entries,  as  described  below. 

Second,  at  most  only  two  initial  conditions  are  required,  wavelength  and  temperature,  and  the  keys 
which  are  used  to  store  these  conditions  may  also  be  use  to  recall  (or  verify)  them.  The  HP41  has  a  flag 
which  it  uses  to  keep  track  of  whether  a  numeric  key  (0  through  9)  has  been  pressed.  This  flag  is  not 
visible  in  the  display.  As  soon  as  any  numeric  key  is  pressed,  this  flag  gets  set.  The  program  uses  this  fact 
to  determine  if  the  user  is  storing  data  or  requesting  an  answer.  If  the  flag  is  set,  the  user  is  assumed  to 
be  storing  data,  and  the  value  in  the  display  is  stored  as  the  new  value  of  the  variable  designated  by  the 
key  pressed.  If  the  flag  is  clear,  the  present  value  of  the  requested  variable  is  displayed.  The  program 
clears  this  flag  automatically  after  every  operation,  so  the  presence  (or  absence)  of  new  numeric  key 
presses  may  be  sensed.  NOTE:  If  this  flag  is  set  and  any  of  the  unshifted  keys  are  pressed,  the  program 
updates  the  wavelength  with  the  current  value  in  the  X  register  (the  display),  in  addition  to  performing 
the  desired  operation.  If  you  have  done  some  simple  keyboard  calculations  (thereby  setting  the  flag),  and 
you  now  wish  to  perform  a  program  operation  without  altering  the  value  for  wavelength,  press  the 
INITIALIZE  key,  first.  This  clears  the  flag  without  altering  the  current  values  for  wavelength  and 
temperature. 

Third,  keys  A,  B,  and  C  are  not  used  for  data  entry,  but  may  be  used  for  configuring  the  program  for 
the  type  of  calculation  to  be  performed.  For  example,  pressing  D  or  E  without  first  pressing  A,  B,  or  C 
will  result  in  the  calculation  of  the  photon  flux  density  in  photons/cmVs/jiM  or  power  density  in 
W/cmVpM,  respectively.  See  the  table  showing  sample  test  data  and  results.  Pressing  A  first,  for  example, 
then  pressing  D  or  E  will  result,  instead,  in  the  calculation  of  the  change  in  flux  (either  photons  or  power) 
per  degree  change  in  blackbody  temperature.  You  may  press  pairs  of  prefix  keys  prior  to  calling  for  a 
calculation  (e.g..  A,  C,  E  to  find  the  change  in  power  up  to  the  cutoff  per  degree  Kelvin).  The  prefix  keys 
may  be  pressed  in  any  sequence  as  long  as  you  end  up  pressing  D  or  E.  NOTE:  It  doesn’t  make  sense  to 
use  both  B  and  C.  The  last  one  pressed  takes  precedence.  If  you  have  pressed  any  prefix  keys  and  change 
your  mind  (or  have  made  an  error),  pressing  the  INITIALIZE  key  will  cancel  and  give  you  a  fresh  start. 
Stored  values  for  wavelength  and  temperature  are  unaffected. 

At  the  end  of  a  calculation,  the  program  displays  in  sequence  (1)  the  key  sequence  which  has  been 
performed,  (2)  the  temperature  that  was  used,  (3)  the  wavelength  that  was  used,  and  (4)  finally,  the 
answer.  The  key  sequence  is  displayed  numerically  proceeding  from  left  to  right,  where  the  numbers 
represent  the  respective  unshifted  keys.  Key  A  is  1,  key  B  is  2;  etc.  For  example,  pressing  C,  A,  D  will 
display  a  key  sequence  of  314.000,  where  the  number  of  trailing  zeros  depends  upon  the  current  display 
setting. 

The  following  equations  and  procedures  are  used  in  BB  il  =  wavelength). 

inf  /  +  p  —  1 

EQN1  =  c,(T/c,)'-p  SUM  ((/ +  p  -  l)!/nJ  EXP(  -  nci/LT)  SUM  „  {(nc2/LT)m/m!)}]. 

n  =  1  m  =  0 

EQN2  =  c. /[L '♦  p ♦ 1  (EXP{c,/LT}  -1)]. 

EQN3  =  c’  EXP  (Cj/LT)  /  IT1  (EXP{  c:/LT}-  1)]. 

EQN4  =  zc,  (T/c,)'  p  (/  +  p-  1)! 


I 
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Default  conditions  are:  p=0,  c,  =  hc7(2Pi)  =  37415,  c;  =  hc/k  =  14388,  zc,  =40495  and  use  EQN2. 

If  prefix  key  A  is  used,  then  p  =  1 ;  if  using  EQN2,  then  multiply  by  EQN3,  and  divide  answer  by  C2. 

If  prefix  key  B  is  used,  then  use  EQN4,  instead. 

If  prefix  key  C  is  used,  then  use  EQN1,  instead. 

If  calculating  power,  then  /  =4. 

If  calculating  photon  flux,  then  /  =3;  use  c,  =  1/k  =  18837E19,  instead;  if  using  EQN4,  then  use 
zCi  =  22642E19,  instead. 

Using  these  equations  and  logical  tests  in  this  fashion  permits  calculation  of  many  different  parameters 
in  a  small  program  space. 


PROGRAM  LISTING  FOR  BB 


01/10/1986 

51  STO  12 

52  RCL  12 

53  RTH 

81*181  *88* 

54*LBL  D 

02*LBL  a 

55  SF  80 

03  37415 

56*LBL  E 

84  STO  04 

57  XEQ  07 

85  14388 

58  RCL  18 

86  STO  10 

59  RCl  11 

07  18837  E19 

68  / 

88  STO  83 

61  STO  13 

89  CIX 

62  RCl  12 

10*LBL  81 

63  X*0? 

11  0 

64  / 

12  STO  15 

65  STO  08 

13  STO  14 

66  4 

14  STO  09 

67  FS?  80 

15  XOY 

68  3 

16  CF  00 

69  STO  85 

17  SF  01 

70  STO  01 

18  CF  22 

71  STO  02 

19  SCI  3 

72  1 

20  RTM 

73  STO  80 

21*LBL  b 

74  F&  81 

22  FS?C  22 

75  GTO  14 

23  STO  11 

76  ST*  01 

24  RCL  11 

77  ST-  02 

25  RTN 

78*181  14 

26H.BL  A 

79  ♦ 

27  CF  01 

88  XEQ  65 

28  XEQ  07 

81  RCl  04 

29  1 

82  FS?  00 

30  GTO  05 

83  RCl  03 

31*181  B 

84  STO  06 

32  SF  02 

85  GTO  I HD  15 

33*LBL  C 

86*LBL  08 

34  XEQ  87 

87  RCl  12 

35  3 

88  RCl  01 

36  FS?C  02 

89  1 

37  2 

90  ♦ 

38  STO  15 

91  YtX 

39*181  05 

92  / 

40  RCL  14 

93  RCl  08 

41  10 

94  EtX 

42  * 

95  STO  89 

43  ♦ 

96  1 

44  STO  14 

97  - 

45  FIX  0 

98  / 

46  RTH 

99  FS?  81 

47*181  07 

108  GTO  01 

48*191  c 

-  101  LASTX 

49  SCI  3 

182  / 

56  FS’C  22 

183  RCL  09 

184  * 

105  RCl  13 

106  Xt2 
10?  * 

188  RCL  18 

109  / 

110  GTO  01 
111*LBL  02 
112  1.202 

113  EHTERt 

114  1.0823 

115  FS?  00 

116  XOY 

117  STO  09 

118  GTO  82 
U9*LBL  03 
128  RCL  01 

121  STO  15 

122  0 

123  RCL  08 

124  RCL  00 

125  * 

126  STO  87 
127*LBL  04 

128  RCL  15 

129  1 

130  - 

131  YtX 

132  LASTX 

133  FACT 

134  / 

135  ♦ 

136  RCL  87 

137  USE  15 

138  GTO  04 

139  EtX 

140  / 

141  RCL  00 

142  RCL  05 

143  YtX 

144  / 

145  ST*  09 

146  RCL  89 

147  / 

148  FS?  08 

149  R-I 

150  FIX  3 

151  RND 

152  ISG  I HD  15 

153  CLfl 

154  X*V 

155  GTO  03 
156*181  82 


157  RCL  06 

158  RCL  13 

159  RCL  02 

168  YtX 
161  / 

162  RCL  01 

163  1 

164  - 

165  FACT 

166  * 

16?  RCl  09 
168  * 

169  FS?  01 

170  GTO  01 

171  RCL  10 

172  / 

173  GTO  01 
I74*LBL  e 

175  .7896 

176  / 
177*LBL  d 

178  3678 

179  / 

180  l/X 

181  SCI  3 

182  CF  22 

183  END 


5.  SV.  This  program  is  a  short,  root-solving  program  which  uses  the  secant  method  to  iterate  to  a  real 
solution  of  arbitrary,  user-defined  functions  of  the  form:  l  (\l  0  It  is  included  here,  since  it  is  used  as  a 
subroutine  by  other  programs  in  this  write-all  configuration.  It  may,  also,  he  used  foi  other  applications 

as  well. 

The  SV  routine  expects  to  find  the  name  of  the  global  label  of  the  function  to  be  solved  in  register  6 
and  an  initial  guess  for  X  in  register  X  (the  display.)  Do  NO  I  use  an  initial  guess  of  zero  for  X.  Registers 

6,  7,  8,  and  9  are  used  by  SV.  XEQing  SV  will  then  result  m  an  iterative,  hopefully  convergent,  search 
for  a  value  of  X  which  makes  the  function  identified  in  register  6  equal  to  zero.  The  successive  values  of 
X  are  displayed,  so  you  can  watch  the  search  converge  (or  diveige).  Lerminatian  of  the  search 
automatically  occurs,  when  two  successive  values  of  \  ate  equal  within  the  resolution  of  the  current 
display  setting.  This  means  that  the  accuracy  of  the  answer  may  be  determined  by  pre-selecting  a  display 
setting.  Upon  termination,  the  solution  for  X  is  in  the  X  register  (display)  and  in  register  7.  NOTE:  When  a 
program  is  running,  the  PRGM  annunciator  is  visible  in  the  display.  When  the  program  stops,  the 
annunciator  goes  out. 

For  example,  suppose  we  wish  to  find  the  solution  to  the  equation: 

X  =  10  M  -  X). 

This  must  be  rewritten  to  take  the  form  F(x)  =  0,  so: 

0  =  X  -  10A(-X). 

A  program  which  performs  this  function  and  which  starts  with  a  global  label  must  be  written.  Any  two- 
letter,  alphabetic  label  is  global,  so  let’s  use  the  label  "FN”.  (Assume  X  starts  in  the  X  register,  and 
the  answer  is  to  be  left  in  the  X  register.)  The  program  would  then  look  like: 

LBL  FN  The  global  label  so  SV  can  call  it. 

ENTER  Starting  value  in  X;  duplicate  it  in  Y. 

10*X  Raise  10  to  the  X  power. 

1/X  Take  the  reciprocal  (yields  10A  (  -  X)). 

Subtract  X  from  Y,  leaving  answer  in  X. 

RTN  Return  to  calling  routine. 

Next,  place  the  label  of  the  function  to  be  solved  into  register  6:  ALPHA  FN  SHIFT  STO  06  ALPHA. 
Set  the  display  to  5  significant  figures:  SHIFT  FIX  5 
Key  in  an  initial  guess:  1 

Execute  the  SV  routine:  XI  Q  Al  PI  LA  SV  ALPHA 

The  routine  will  be  observed  to  go  through  six  iterations;  finally,  stopping  at  a  solution  for  X  of 
.39901.  Executing  FN  with  this  value  in  the  X  register  produces  a  result  of  1  E-10  which  is  reasonably 

close  to  zero. 

SV  properly  converges  for  routines  in  this  WAIT  .  However,  lor  some  user  functions,  it  may  have 
problems  with  bad  initial  guesses,  or  where  multiple  solutions  exist  (particularly,  when  some  are 
imaginary),  or  when  Xs  near  a  root  are  invalid  arguments  (e  g.SORt.X)  0  can't  tolerate  negative  ,\s.) 


PROGRAM  LISTING  FOR  SV 


08  OUG  86 


PR  P 

••SV" 

01 ♦LBL 

••  s  v  •• 

82 

STO 

07 

03 

XEQ 

IND 

06 

04 

STO 

08 

05 

RCL 

07 

06 

1 

07 

08 

STO 

09 

09 

ST- 

07 

10*LBL 

00 

1  1 

VIEW  07 

12 

RCL 

07 

13 

XEQ 

IND 

06 

14 

ENTERf 

15 

ST- 

08 

16 

X<  > 

08 

17 

/ 

18 

ST* 

09 

19 

RCL 

09 

20 

ST- 

07 

21 

RCL 

07 

22 

+ 

23 

RND 

24 

RCL 

07 

25 

RND 

26 

X*Y? 

27 

GTO 

00 

28 

RCL 

07 

29 

RTN 
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6.  GR-DIFF. 


SHIFT  fl  SHIFT  B  SHIFT  C  SHIFT  D  SHIFT  E 


ENTER 

or 

FIND 
wavelength 
for  present 
condition 


ENTER 

ENTER 

INTERCHANGE 

TOGGLE 

or 

or 

CONDITION 

GR/DIFF 

FIND 

FIND 

1  AND  2 

detector  temp 

R  or  RoA 

for  present 

for  present 

condition 

condition 

Below  is  a  table  showing  six  example  pieces  of  data  describing  two  different  test  conditions.  Any  five  of 
these  may  be  entered,  and  the  program  will  solve  for  the  sixth.  The  solution  will  assume  either  a  diffusion 
limited  detector  or  a  G-R  limited  detector,  depending  on  the  status  of  the  GR-DIFF  toggle. 


CONDITION 

WAVELENGTH  (pM) 

TEMPERATURE  (K) 

R  or  RoA 

1 

9 

80 

10 

2  (DIFF  LTD.) 

11.5 

65 

4.8 

2  (G-R  LTD.) 

11.5 

65 

6.9 

NOTE:  If  you  can’t  remember  the  key  definitions,  press  SHIFT  RTN  R/S. 

The  purpose  of  this  program  is  to  permit  you  to  enter  reported  performance  data  for  detectors  which 
have  the  “wrong”  cutoff  or  were  tested  at  the  “wrong”  temperature  and  to  extrapolate  how  the  detectors 
would  perform,  if  they  had  the  “right”  cutoff  and/or  were  tested  at  the  “right”  temperature.  Or,  for 
example,  you  can  take  reported  data  and  predict  what  temperature  and/or  wavelength  of  operation  would 
be  required  to  meet  a  desired  R  or  RoA.  One  must  assume  a  measure  of  quality  of  the  detectors  by 
specifying  whether  they  are  diffusion  limited  or  G-R  limited,  this  is  done  by  pressing  the  GR/DIFF 
toggle.  Current  status  is  indicated  prior  to  any  result  being  displayed. 


Only  three  keys  (A,  B,  and  C)  are  used  to  enter  any  of  the  six  data.  A  is  used  for  wavelength  for  both 
test  conditions.  B  is  used  for  temperature  for  both  test  conditions.  C  is  used  for  R  or  RoA  for  both  test 
conditions.  Alternating  between  conditions  is  done  with  the  D  key.  After  any  entry,  the  display  will  show 
the  status  of  ALL  THREE  parameters  for  the  current  test  condition.  Read  these  values  from  left  to  right 
in  the  same  order  as  keys  A,  B,  and  C.  That  is,  wavelength,  temperature,  R  or  RoA.  The  parameters  are 
separated  by  a  ?  mark  until  a  calculation  is  requested  which  makes  all  six  data  consistent. 

For  example,  viewing  the  sample  data  at  the  beginning  of  par  6.,  suppose  we  were  presented  with  test 
condition  1  and  wanted  to  know  what  RoA  to  expect  for  a  diffusion  limited  diode  under  test  condition  2. 
If  you’re  not  sure  whether  the  GR-DIFF  toggle  is  set  to  the  diffusion  limited  case,  press  the  toggle  once 
and  see  what  condition  occurs.  If  it’s  OR;  then,  press  it  again.  NOTE:  Performing  a  CLRG  (SHIFT 
TAN)  clears  all  registers  and  defaults  to  the  diffusion  limited  case.  Press  9,  A,  80,  B,  10,  C,  D  (to 
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interchange  which  condition  keys  A,  B,  and  C  affect,  and  which  condition  is  displayed),  11.5,  A,  65,  B, 

C  (to  solve  for  RoA  for  this  condition).  You  may  toggle  between  (and  view)  the  two  test  conditions,  at 
any  time,  by  pressing  D.  You  may  toggle  between  diffusion  limited  and  G-R  limited  operation,  at  any 
time,  by  pressing  E.  NOTE:  If  the  calculation  has  been  performed  for  one  mode  of  operation  (Say,  G-R), 
and  the  other  mode  is  selected,  ?  marks  are  restored  to  the  display  to  indicate  inconsistent  data. 

The  program  may  also  be  used  when  one  of  the  test  conditions  is  diffusion  limited,  and  the  other  is 
G-R  limited.  But  you  have  to  get  from  here  to  there  in  two  steps,  and  you  must  know  (or  assume)  the 
temperature  at  which  the  break  from  diffusion  limited  to  G-R  limited  operation  occurs.  For  example, 
suppose  a  diode  with  a  4.2  pM  cutoff  at  195  K  has  an  Ro  of  1E6  ohms.  The  break  from  diffusion  limited 
to  G-R  limited  operation  for  these  diodes  normally  occurs  near  130  K.  You  want  to  know  what  Ro  to 
expect  at  80  K.  Strictly  speaking,  the  cutoff  will  have  changed  when  the  diode  is  cooled  to  80  K,  but 
ignore  that  fact  here.  The  procedure  is  as  follows:  CLRG  (clears  all  registers  and  puts  the  machine  into 
diffusion-limited  regime,  4.2,  A,  195,  B,  1E6,  C  (test  condition  1  now  contains  the  195  K  performance), 
D,  4.2,  A,  130,  B,  C  (test  condition  2  now  contains  the  Ro  at  130  K,  diffusion  limited,  Ro  =  6.5  E9),  E 
(to  toggle  to  G-R  limited),  D  (test  condition  1  will  now  represent  the  80  K  performance)  (wavelength  is 
still  4.2,  so  no  need  to  reenter),  80,  B,  C  (to  Find  the  G-R  limited  Ro  at  80  K).  The  answer  is  Ro  =  2.5 
E13.  More  correct  would  be  to  use  the  correct  predicted  values  for  the  wavelength  at  130  K  and  80  K,  but 
solution  for  such  parameters  won’t  be  covered  until  discussion  of  program  HCT  (par  7.). 


PROGRAM  I  I  s  M  M «  FOR  GR-DIFF 


01/18/1986 


81*LBL  “GR— BIFF” 

92  •«, T.Rfl. SHIFT, GS- 
03  PROMPT 
04*LBL  D 
85  RCL  89 
06  X<>  10 
07  STO  09 
98  RCL  08 
09  X<>  11 

10  STO  98 

11  RCL  07 

12  X<>  12 

13  STO  07 
14*LBL  81 

15  CL  A 

16  FIX  2 

17  ARCL  99 

18  ARCL  14 

19  FIX  8 

20  flRCL  08 

21  flRCL  14 

22  SCI  1 

23  flRCL  97 

24  PROMPT 
25*LBL  04 
26  1 

27  RCL  13 

28  X=Y? 

29  RTN 
38  2 

31  X=Y? 

32  P  TH 
33*LBL  E 

34  •?' 

35  flSTO  14 

36  RCL  13 

37  1 

38  X=Y? 

39  2 

48  STO  13 
41 ♦LSL  02 

42  -DIFFUS.H- 

43  RCL  13 

44  2 

45  X=Y? 

46  -G-R* 

47  -F  LTD- 
43  fiV  TEH 

49  RTH 
50*LBL  C 


51  ->• 

52  flSTO  14 

53  STO  07 

54  FS?C  22 

55  GTO  81 

56  XES  84 

57  XES  02 
53  14388 

59  RCL  13 

60  / 

61  STO  86 

62  RCL  99 

63  RCL  08 

64  * 

65  / 

66  RCL  06 

67  RCL  18 

68  RCL  11 
59  * 

78  / 

71  - 

72  EtX 

73  RCL  12 

74  * 

75  STO  87 

76  -  * 

77  flSTO  14 

78  GTO  81 
79»LBL  B 
80  •?* 

81  flSTO  14 

82  STO  08 

83  FS?C  22 

84  GTO  81 
35  GTO  08 
86*LBL  A 
87  *?' 

83  flSTO  14 

89  STO  89 

90  FS?C  22 

91  GTO  01 

92  SF  80 
93»LBL  08 

94  XEO  04 

95  XES  02 
95  RCL  87 

97  PCL  12 

98  / 

99  LH 

100  14338 

101  / 

102  RCL  13 

1  07  t 


184  RCL  18 

105  RCL  11 

106  * 

187  1/X 
108  ♦ 

109  1/X 

110  RCL  09 

111  RCL  88 

112  FS?  80 

113  XOY 

114  RUH 

115  / 

116  FS?  80 

117  STO  09 
113  FC?C  00 

119  STO  08 

120  ‘  ‘ 

121  flSTO  14 

122  CTO  01 

123  EHB 


7.  HCT,  X,  EG. 


Registers  used:  4  —  Cutoff  wavelength  OR  bandgap 

5  —  Temperature 

6  —  “X”  alpha  label  for  SV  program 

7  —  X,  either  from  direct  entry  or  left  by  SV  program 

8  —  Used  by  SV 

9  —  Used  by  SV 

This  program  uses  the  empirical  equation  for  bandgap  in  terms  of  temperature  and  X  value  developed 
by  Hansen,  Schmidt,  and  Casselman  (HSC)',  which  is  valid  for  a  wide  range  of  X: 

Eg  =  -  .302  +  1.93  X  +  (5.35E-4)T(1  -  2X)  -  0.810  X  A2  +  0.832  XA  3. 

It  is  generally  desirable  to  determine  the  bandgap  or  cutoff  that  a  device  would  have  at  some 
temperature  from  data  which  are  available  at  some  other  temperature.  The  parameter  which  remains  fixed 
in  such  calculations  is  the  material  composition,  X.  If  X  is  known,  it  is  stored  into  the  program  via  label 
C;  labels  D  and  E  calculate  bandgap  or  cutoff  directly  from  temperature.  Sometimes,  however,  the  data 
provided  does  not  include  X,  so  the  HSC  equation  is  solved  iteratively  for  X  as  a  function  of  T  and  (Eg  or 
cutoff)  through  the  use  of  labels  A  or  B.  The  value  found  for  X  is  displayed  and  saved  by  the  program 
for  use  by  D  and  E. 

This  program  solves  the  HSC  equation  for  X  by  calling  the  SV  program,  described  elsewhere  (par  5.), 
and  searching  for  a  root  of  the  equation  HSC  -  Eg  =  0. 

Global  labels  available  from  this  program  are  HCT  which  behaves  like  local  label  A  and  may  be  used 
to  find  X  from  cutoff  and  temperature,  and  EG  which  calculates  Eg  from  an  internally  stored  value  for  T 
and  an  X  in  the  X  register.  Of  course,  these  may  be  called  from  other  programs. 

'Hansen,  Schmidt,  and  Casselman  (HSC),  J.  Appl.  Phys.  53,  7099;  1982. 
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PROGRAM  LISTING  FOR  HCT.  X,  EG 


88  HUG  86 


81*LBL  'HCT 
02»L8L  ft 
63  CF  80 
04  CTO  01 
85*L8L  B 
06  SF  08 
87*L8L  81 
08  FIX  9 
09  *X* 

10  ASTO  06 

11  STO  05 

12  RON 

13  STO  04 

14  CLX 

15  ENTERt 

16  .25 

17  SF  10 

18  XEQ  *SV* 

19  RTN 
28*LBL  E 

21  XEfl  D 

22  1.242 

23  / 

24  1/X 

25  RTN 
264LBL  -X* 

27  XEQ  00 

28  FC?  00 

29  1.242 
38  RCL  04 

31  FC?  80 

32  / 

33  - 


34  RTN 
35*LBL  D 

36  STO  05 

37  RCL  07 
38*LBL  06 
39*LBL  "EG- 

40  ENTERt 

41  ENTERt 

42  ENTERt 

43  .832 

44  * 

45  .81 

46  - 

47  * 

48  1.93 

49  + 

50  5.35  E-4 

51  RCL  85 

52  * 

53  ENTERt 

54  RON 

55  ST+  X 

56  - 

57  * 

58  ♦ 

59  .382 

60  - 
61  RTN 
62*LBL  C 

63  STO  07 

64  END 
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8.  NI,  R0A. 


majority 


lifetime 

carrier 

X 

cutoff 

detector 

(sec) 

concentration 

(/cc) 

(composition) 

(mM) 

temperature 

(K) 

SHIFT  R 

SHIFT  B 

SHIFT  C 

SHIFT  D 

SHIFT  E 

base  layer 

detector  area 

FIND 

FIND 

FIND 

thickness 

(sq.  cm.) 

saturation 

RoA  product 

intrinsic 

0*M) 

current 

(ohm*sq.cm.) 

carrier 

(amps)  concentration 

(/sq.cm.) 


Below  is  a  table  showing  the  data  required,  registers  used,  sample  data,  and  results. 


MINORITY  CARRIER  LIFETIME  (s) 

MAJORITY  CARRIER  CONCENTRATION  (/cm!) 
COMPOSITION  (X)  OR 

CUTOFF  WAVELENGTH  O'M) 

DETECTOR  TEMPERATURE  (K) 

BASE  LAYER  THICKNESS  0<M) 

DETECTOR  AREA  (cm*)  (only  req’d  for  Isat) 

TEST  DATA 

IE-7 

5E15 

.21200 

12 

80 

10 

6.45E-6 

REGISTER  USED 

1 

2 

7 

4 

5 

3 

12 

RESULTS  FOR  CUTOFF  =  12 

INTRINSIC  CARRIER  CONCENTRATION  (/cm1) 

3.55 1 E 1 3 

PREDICTED  RoA  (ohm'cm1) 

17.064 

SATURATION  CURRENT  (A) 

2.606E-9 

This  program  uses  the  following  equation  to  solve  for  RoA.  It  assumes  that  the  device  is  diffusion 
limited,  that  the  diffusion  length  is  much  longer  than  the  base  layer  is  thick,  and  that  surface  effects  are 
not  present.  Also,  lateral  collection  effects  are  not  considered,  unless  one  takes  the  A  as  the  effective 
carrier  collecting  area  and  not  the  junction  area. 

RoA  =  kTNat  , 
e*2  d  NiA2 

where  e  is  the  charge  on  an  electron; 
t  is  the  minority  carrier  lifetime; 
k  is  Boltzmann’s  constant; 

T  is  the  detector  temperature; 
d  is  the  base  layer  thickness; 

Na  is  the  majority  carrier  concentration;  and 
Ni  is  the  intrinsic  carrier  concentration. 
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All  these  parameters  are  user  input,  except  the  intrinsic  carrier  concentration  which  must  be  solved  for 
from  the  detector  temperature  (already  input),  and  the  material  composition  (X).  The  equation  used  is 
from  Hansen  and  Schmidt2: 

Ni  =  [5.585  - 3.820 X  +  (1 .753E-3)T  - (1 ,364E-3)XT]  [  (10“)Eg75T'  VE8/2kT). 

You  have  the  option  of  providing  the  materia!  cutoff  wavelength  instead  of  X,  in  which  case,  the 
program  calls  the  HCT  program  to  solve  for  X,  and  uses  (1.24/cutoff)  for  Eg.  If  X  is  provided  instead  of 
the  cutoff,  the  program  calls  the  EG  program  to  find  Eg.  NI  keeps  track  of  whether  X  or  the  cutoff  has 
been  input  (and  selected  to  predominate.)  For  example,  if  calculations  are  done  for  various  device 
temperatures  and  X  predominates,  new  values  for  Eg  will  be  calculated  at  each  temperature.  If  cutoff 
predominates,  new  values  of  X  will  be  calculated  at  each  temperature.  If  it  is  desired  that  the  X  value 
predominate,  but  only  the  cutoff  is  known,  run  the  program  once  entering  the  known  cutoff.  Then, 
EITHER  quickly  jot  down  the  X  value  that  is  briefly  displayed  and  key  it  in  via  SHIFT  C  prior  to  any 
further  calculations,  OR  merely  clear  flag  2. 

The  diode  diffusion  limited  saturation  current  may  be  found,  once  the  RoA  is  known,  IF  the  diode 
area  is  known.  The  equation  used  is: 


Isat  =  k  T  ,  where  Ro  =  Ro  A. 
e  Ro  A 

As  a  point  of  interest,  viewing  the  program  code  for  LBL  C  (Isat),  one  observes  that,  despite  the  fact 
that  R0A  should  be  executed  as  a  subroutine,  the  ROA  routine  is  not  run  via  an  XEQ  (which  is  a 
subroutine  call)  but  via  a  GTO  (which  is  a  jump);  and  that  ROA  returns  to  the  Isat  routine  based  on  the 
status  of  flag  3.  The  reason  for  this  more-involved  method  for  calling  ROA  lies  in  the  fact  that  the 
HP41’s  RTN  (return)  stack  is,  only,  six  levels  deep.  This  stack  is  where  the  HP41  keeps  track  of  the 
location  from  which  a  subroutine  has  been  called  (the  location  where  it  must  return  at  the  end  of  a 
subroutine.)  Following  the  program  code,  one  finds  that  ROA  calls  NI,  which  may  call  label  01,  which 
calls  HCT,  which  calls  SV,  which  calls  X,  which  calls  EG.  If  one  tries  to  use  ROA  as  a  subroutine,  the 
return  address  in  the  initial  calling  program  gets  pushed  off  the  stack  (and  lost),  when  EG  gets  called;  so 
the  program  terminates  (incorrectly)  in  ROA. 


2Hansen  and  Schmidt,  J.  Appl.  Phys.  54,  3;  March  1983. 


PROGRAM  LISTING  FOR  NI,  R0A 


01/10/1984 

51*LBL  92 

183*LBL  a 

52  RCL  07 

184  STO  81 

01*181  "HI* 

53  XEQ  'EG* 

105  RTH 

02*181  E 

54  1.24 

106*LBL  b 

03  RCL  11 

55  / 

107  STO  82 

04  FS?  01 

56  1/X 

108  RTH 

65  RTN 

57  FIX  2 

109*LBL  c 

04  FS?  82 

58  *<i=* 

118  STO  87 

07  XEQ  01 

59  fiRCL  X 

111  CF  01 

88  FC?  02 

60  AVIEH 

112  CF  82 

09  XEQ  02 

61  STO  04 

113  RTH 

10  -1.364 

62  RCL  07 

114*LBL  d 

11  RCL  05 

63  RTH 

115  STO  84 

12  * 

64*LBL  ■R0&* 

116  CF  81 

13  3820 

65*LBL  D 

117  SF  82 

14  - 

66  FIX  3 

118  RTH 

15  * 

67  XEQ  E 

119*LBL  e 

16  1.753 

68  5378  E15 

128  STO  85 

17  RCL  05 

69  XOY 

121  CF  81 

18  * 

78  Xt2 

122  RTH 

19  + 

71  / 

123*LBL  A 

20  5585 

72  RCL  05 

124  STO  03 

21  * 

73  * 

125  END 

22  Ell 

74  RCL  02 

23  * 

75  * 

24  1.24 

76  RCL  03 

23  RCL  04 

77  / 

26  / 

78  RCL  01 

27  SORT 

79  • 

28  RCL  85 

88  -Rfl=* 

29  * 

81  ARCL  X 

38  SORT 

82  AVIEH 

31  3 

83  FS?C  03 

32  YtX 

84  GTO  03 

33  * 

85  RTH 

34  7194 

86*LBL  8 

35  RCL  04 

87  STO  12 

36  / 

88  RTN 

37  RCL  85 

89*LBL  C 

38  / 

90  SF  83 

39  EtX 

91  GTO  D 

40  / 

92*LBL  03 

41  SF  01 

93  RCL  12 

42  STO  11 

94  / 

43  CLD 

95  RCL  05 

44  FrX  3 

96  / 

45  RTH 

97  11605 

46*L8L  01 

98  * 

47  RCL  04 

99  1/X 

48  RCL  05 

100  -I*’ 

49  XEQ  *HCT* 

101  QRCL  X 

50  RTH 

102  PROMPT 

9.  CONSTNT. 


c  eo  sigma  R  N* 


speed  of 
light 

permitivity 
of  free 
space 

Stefan- 

Boltzmann 

constant 

gas 

constant 

Avogadro’s 

number 

SHIFT  R 

SHIFT  B 

SHIFT  C 

SHIFT  D 

SHIFT  E 

k 

Boltzmann 

constant 

e 

charge  on  an 
electron 

h 

Plank 

constant 

h/(2*PI) 

Plank 

constant 

mt 

mass  of  an 
electron 

divided  by 
(2*PI) 


This  program  may  be  used  to  conveniently  bring  any  of  the  available  physical  constants  into  the  X 
register  for  use  in  direct  keyboard  calculations.  All  constants  are  in  MKS  units.  When  any  of  the  keys  are 
pressed  (SHIFT  A  through  SHIFT  E,  or  A  through  E)  the  current  value  in  the  X  register  is  pushed  up 
into  the  Y  register,  and  the  selected  constant  is  placed  in  the  X  register.  For  example,  suppose  one  wants 
to  evaluate  the  following  expression: 

h  c  ,  where  wavelength  =  4.2E-6  m, 

wavelength  k  T  and  T  =  300  K. 

The  sequence  that  may  be  followed  is: 


PRESS 

TO  OBTAIN 

SEE 

SHIFT  SCI  4 

Floating  point  display  and  5  significant  figures 

C 

h 

6.6262  -34 

SHIFT  A 

c 

2.9979  08 

* 

Multiply  X  and  Y  registers 

1.9865  -25 

4.2  EEX  CHS  6 

Wavelength 

4.2  -6 

/ 

Divide  Y  register  by  X 

4.7297  -20 

A 

k 

1.3806  -23 

/ 

Divide  Y  register  by  X 

3.4258  03 

300 

T 

300 

/ 

Divide  Y  register  by  X 

1.1419  01 
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PROGRAM  LISTING  FOR  CONSTNT 


01/18/1986 

01*IBL  -COHSTNT* 
02*LBL  a 
03  2.99792S  E8 
04  RTH 
0SKBL  fl 
06  1.38062  E-23 
07  RTH 
88*IBL  B 
09  1.6022  E-19 
10  RTH 
11M.BL  C 

12  6.6262  E-34 

13  RTH 
14*UJL  0 

15  1.0546  E-34 

16  RTR 
17*LBL  E 

18  9.1896  E-31 

19  RTH 
20*LBL  e 

21  6.02217  E23 

22  RTH 
23*LBL  b 

24  8.85  E-12 

25  RTN 
26*LBl  c 

27  5.66961  E-8 

28  RTH 
29H.BL  d 
38  8.31434 


10.  IN' JEFF. 


detector 

area 

(sq.  cm.) 

wavelength 
cuton  OiM) 
ENTER 
cutoff  (>iM) 

background 

temperature 

(K) 

quantum 

efficiency 

(percent) 

f  nt 

(cold  shield) 

SHIFT  H 

SHIFT  B 

SHIFT  C 

SHIFT  D 

SHIFT  E 

zero  bias 
resistance 
(ohms) 


detector  bias 
(volts) 


detector 

temperature 

(K) 


detector  m 
value  (1  = 
diff.  ltd) 
ENTER 
CCD  input 
derating 
factor 

(usually  =  1.4) 

D 


FIND 

Rd 

Q 

lp 

Iv 

Rf 

IE 

E 


Below  is  a  table  showing  the  data  required,  registers  used,  sample  data,  and  results: 


PARAMETER 


REGISTER  USED  SAMPLE  DATA  RESULTS 


DETECTOR  AREA,  A 

16 

IE-5 

CUTON 

5 

3.0 

CUTOFF 

10 

5.0 

BACKGROUND  TEMPERATURE 

12. 

300 

QUANTUM  EFFICIENCY 

17 

80 

COLD  SHIELD  F/tt 

13 

1 

ZERO  BIAS  RESISTANCE,  Ro 

9 

1E6 

DETECTOR  BIAS,  V 

2 

-  .02 

DETECTOR  TEMPERATURE 

14 

200 

DETECTOR  m  VALUE  (1  =  DIFF.  LTD.,  2  =  G-R  LTD.) 

3 

1.2 

CCD  DERATING  FACTOR  (m  FOR  CCD  INPUT  NODE) 

11 

1.4 

Rd  =  DIODE  DYNAMIC  RESISTANCE  (OHMS) 

1 

2.63E6 

Q  =  PHOTON  FLUX  (/cmVs) 

15 

1.31E16 

Ip  =  DIODE  PHOTOCURRENT  (A) 

4.20E-9 

Iv  =  DIODE  BIAS  CURRENT  (A) 

I.28E-8 

Rf  =  CCD  INPUT  IMPEDANCE  (OHMS) 

1 .42F.6 

(RECIPROCAL  OF  TRANSCONDUCTANCE) 

IE  =  INJECTION  EFFICIENCY 

0.650 

Other  registers  used:  4,  6,  7,  8 


This  program  wiiculaie^  injection  efficiency  from  the  equation: 


IE  =  _ Rd _ . 

(CCD  derating  factor)Rf  +  Rd 

Rf  is  it'.*  c  •;  D  input  impedance  and  is  the  reciprocal  of  the  term  normally  used:  CCD  input 
transconductance.  Rf  is  computed  from: 

Rf  =  k  T  ,  where  Id  is  the  net  diode  current  =  Iv  +  Ip. 

Tld 

Ip  is  the  photocurrent  induced  in  the  photodiode  by  the  background  photon  flux,  Q,  and  is: 

Ip  =  e  A  Q  (quantum  efficiency). 

Iv  is  the  diode  current  due  to  the  fact  that  a  bias  is  applied  and  is  found  from  the  equation  of  a  diode: 

Iv  =  mkT  [EXP(  -  eV/{mkT))-  1). 
e  Ro 

Finally,  Rd  is  the  dynamic  resistance  of  the  diode  at  the  specified  bias.  This  may  be  found  by 
differentiating  the  equation  for  Iv  with  respect  to  V  and  taking  the  reciprocal,  so  that: 

Rd  =  RoEXP(-eV/{mkT}). 

When  E  is  pressed  to  compute  an  injection  efficiency,  the  intermediate  results  (Rd,  Q,  Ip,  Iv,  Rf),  with 
identification,  are  flashed  in  the  display,  as  they  are  found.  The  program  stops  with  the  injection 
efficiency  displayed.  If  a  printer  is  connected  to  the  HP  41  and  is  on  before  the  HP41  is  turned  on,  all 
results  will  be  printed.  These  data  flash  past  too  fast  to  be  recorded  without  a  printer  and,  if  you’re  not 
interested  in  it,  just  ignore  it.  If  your  ARE  interested  in  the  intermediate  results,  set  flag  21  using  SHIFT 
SF  21  prior  to  pressing  E.  The  HP  41  will  stop  at  each  intermediate  value  so  that  it  may  be  recorded. 
Pressing  R/S  at  each  of  these  points  will  continue  the  program  to  the  next  result,  ultimately  displaying  the 
value  for  the  injection  efficiency.  You  need  set  flag  21  only  once,  and  it  will  remain  set  until  you  clear  it 
using  the  FLAG  routine,  using  SHIFT  CF  21,  or  by  turning  the  HP4I  on  with  no  printer  connected. 


I 

i 
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PROGRAM  LISTING  FOR  INJEFF 


01/10/1986 

51  "Rd=* 

184  RCL  01 

52  ARCL  X 

185  ♦ 

53  AVIEH 

186  RCL  81 

01*LBL  MNJEFf 

54  FS’  90 

107  / 

02*LBl  a 

55  GTO  00 

108  1/X 

03  STO  16 

56  CF  02 

109  *IE=* 

04  RTH 

57  RCL  10 

118  FIX  3 

05*LBL  b 

58  XEQ  -QBB- 

111  ARCL  X 

06  STO  10 

59  STO  15 

112  AVIEH 

07  XOY 

60  RCL  05 

113  EK9 

08  STO  85 

61  XEQ  -Q8B- 

09  CF  00 

62  ST-  15 

10  RTH 

63  SF  90 

ll*LBl  c 

64*IBL  00 

12  STO  12 

65  RCL  15 

13  CF  00 

66  -8=* 

14  RTH 

67  ARCL  X 

15*LBl  d 

68  AVIEH 

16  STO  17 

69  RCL  17 

17  RTH 

70  * 

18*LBL  e 

71  RCL  13 

19  STO  13 

72  Xt2 

20  RTH 

73  / 

21*LBL  A 

74  4  E-22 

22  STO  09 

75  * 

23  RTH 

76  RCL  16 

24*181  B 

77  ♦ 

25  STO  02 

78  CHS 

26  RTH 

79  *IP=' 

27*LBL  C 

88  ARCL  X 

28  STO  14 

81  AVIEH 

29  RTH 

82  RCL  84 

38HJ5L  D 

83  1 

31  STO  11 

84  - 

32  XOY 

85  RCL  06 

33  STO  03 

86  / 

34  RTH 

87  RCl  89 

35*LBL  E 

88  / 

36  SCI  2 

89  RCL  03 

37  11605 

90  • 

38  RCl  14 

91  *1V=* 

39  / 

92  ARCL  X 

40  STO  06 

93  AVIEH 

41  RCl  02 

94  ♦ 

42  * 

95  ABS 

<3  RCl  03 

96  RCL  06 

44  / 

97  * 

45  EtX 

98  1/X 

46  STO  84 

99  RCL  11 

47  t/X 

108  * 

48  RCl  09 

181  -RF=' 

49  * 

102  ARCL  X 

50  STO  01 

103  AVIEH 

30 

III.  FINAL  COMMENTS 


As  with  any  software,  the  results  are  only  as  valid  as  the  data,  the  equations,  and  the  program  steps 
which  are  intended  to  represent  them.  The  various  programs  have  been  tested  and  are  believed  to  be 
correct  and  functional.  However,  if  any  errors  are  found,  if  any  problems  arise,  or  if  you  have  and 
comments  or  questions,  please  let  the  author  of  this  report  know. 
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1  Commander 
ARRADCOM 
ATTN:  DRDAR-SEC 
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1  Commander 
CORADCOM 
ATTN:  DRDCO-PPA-CA 
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1  Commander 
CERCOM 
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1  Commander 
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1  Director 
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Alexandria,  VA  22314 

1  Project  Manager 
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1  Project  Manager 
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1  Commandant 

US  Army  Air  Defense  School 
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US  Army  Armor  School 
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of 
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of 

Copies 


Addressee 


1  Honeywell,  Inc. 
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